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ABSTRACT 

The progenitor of type Il-plateau supernova (SN) 2004dj is identified with a supergiant in a compact 
star cluster known as "Sandage Star 96" (S96) in the nearby spiral galaxy NGC 2403, which was 
fortuitously imaged as part of the Beijing- Arizona- Taiwan-Connecticut (BATC) Multicolor Sky Survey 
from Feb 1995 to Dec 2003 prior to SN 2004dj. The superior photometry of BATC images for S96, 
taken with 14 intermediate-band filters covering 3000 — 10000A, unambiguously establishes the star 
cluster nature of S96 with an age of ~ 20Myr, a reddening of Fi(B — V) ~ 0.35 mag and a total mass of 
~ 96, OOOMq. The compact star cluster nature of S96 is also consistent with the lack of light variations 
in the past decade. The SN progenitor is estimated to have a main-sequence mass of ~12Mq. The 
comparison of our intermediate-band data of S96 with the post-outburst photometry obtained as the 
SN has significantly dimmed, may hopefully conclusively establish the nature of the progenitor. 
Subject headings: galaxies: individual (NGC 2403) — galaxies: star clusters — stars: evolution — 
supergiants — supernovae: general — supernovae: individual (SN 2004dj) 
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1. INTRODUCTION 

Identification of the progenitors of SNe is extremely 
valuable for testing theories of stellar evolution and 
supernova explosions. Type II SNe arise from core- 
collapses of evolved massive stars (Paczynski 1971; Gol- 
dreich & Weber 1980). The subclass of type II-P SNe is 
thought to be associated with red supergiants of higher 
initial masses that retained their hydrogen envelopes be- 
fore core collapse. This model accounts for the main 
observed features (spectra and light curves) and the esti- 
mated physical parameters of the expanding photosphere 
such as velocity, temperature and density (Woosley & 
Weaver 1986; Hamuy 2003). However, there is little di- 
rect evidence for the red supergiant hypothesis. To date, 
only three SNe have had their progenitors identified: SN 
1987A (Gilmozzi et al. 1987; Sonneborn et al. 1987), SN 
1993 J (Aldering et al. 1994; Maund et al. 2004), and 
SN 2003gd (Van Dyk et al. 2003; Smartt et al. 2004). 
The progenitor of peculiar type II SN 1987A was a blue 
supergiant. The Type lib SN 1993J arose in a massive in- 
teracting binary. The progenitor of type II-P SN 2004et 
in NGC 6946 was tentatively identified as a yellow super- 
giant (Li et al. 2005). Thus the expected red supergiant 
origin for the common type II-P SNe is only favored for 
SN 2003gd. The fortuitous occurrence of type II-P SN 
2004dj in a nearby galaxy with prediscovery images has 
allowed us to examine its progenitor. 

SN 2004dj was discovered on July 31.76 UT 2004 by K. 
Itagaki (Nakano 2004) in the nearby spiral galaxy NGC 
2403 about 3.3 Mpc away (Karachentsev et al. 2004). 
Initially, the reported V— band magnitude was ~ 11.2 
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mag, making SN 2004dj the optically brightest SN in the 
past decade since SN 1993J in M81. The initial spec- 
trum of SN 2004dj resembles that of a normal type II-P 
supernova, with prominent P-Cygni profiles in hydrogen 
Balmer lines (Patat et al. 2004). 

Through astrometric registration of archival images 
of NGC 2403 and recent images of SN 2004dj, Mai'z- 
Apellaniz et al. (2004; MA04) established that SN 2004dj 
coincided with Sandage Star 96 (S96) in the list of lumi- 
nous stars and clusters in NGC 2403 (Sandage 1984). Us- 
ing the earlier photometry of S96 (Larsen 1999), MA04 
suggested that S96 was a young compact cluster with an 
age of 13.6 Myr and a total stellar mass of 24,000 Mq. 
They inferred that the progenitor of SN 2004dj had a 
main-sequence mass of ~ 15M Q . 

In this Letter, we identify the progenitor of SN 
2004dj by combining recent and archival images from the 
Beijing-Arizona- Taiwan-Connecticut (BATC) Sky Sur- 
vey. We find that the spectral energy distribution (SED) 
of S96 resembles that of a star cluster rather than a single 
star. Using the simple stellar population (SSP) model, 
we re-examine the age and mass of S96 and set new limits 
on the progenitor mass of SN 2004dj . 

2. BATC OBSERVATIONS OF SANDAGE STAR 96 

2.1. Archival Images of the BATC Sky Survey 

The observations of NGC 2403 were obtained by the 
BATC 60/90cm Schmidt telescope located at the Xing- 
Long station of the National Astronomical Observatory 
of China (NAOC). This telescope has 15 intermediate- 
band filters covering the optical wavelength range of 
3000 — lOOOOA, and is specifically designed to avoid con- 
taminations from the brightest and most variable night- 
sky emission lines. Descriptions of the BATC photomet- 
ric system can be found in Fan et al. (1996). 

Figure 1 compares the pre- and post-explosion images 
of SN 2004dj. The left panel shows the field of SN 2004dj 
and neighboring bright field stars, as imaged in the i 
band (centered on 6656 A) of the BATC system with 
the NAOC 60/90cm Schmidt telescope on 2004 Aug 11. 
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TABLE 1 

BATC Photometry of Sandage Star 96 



Filter A(A) FWHM(A) 


N a 


Value(mag) 


Time Span 


a 


3360 


360 


8 


18.73(0.24) 


Dec 1998 


c 


4210 


320 


8 


18.18(0.04) 


Dec 1998 - 


Nov 2003 


d 


4540 


340 


38 


18.16(0.03) 


Dec 1998 - 


Nov 2003 


e 


4925 


390 


17 


18.12(0.07) 


Mar 1995 - 


Jan 2003 


f 


5270 


340 


9 


18.12(0.05) 


Feb 1995 - 


Dec 1995 


g 


5795 


310 


13 


17.93(0.04) 


Feb 1995 - 


Jan 1996 


h 


6075 


310 


12 


17.85(0.06) 


Feb 1995 - 


Jan 1996 


i 


6656 


480 


11 


17.78(0.06) 


Feb 1995 - 


Jan 1999 


j 


7057 


300 


10 


17.61(0.09) 


Feb 1995 - 


Dec 2003 


k 


7546 


330 


13 


17.57(0.06) 


Feb 1995 - 


Mar 1996 


m 


8023 


260 


17 


17.57(0.08) 


Feb 1995 - 


Dec 2003 


n 


8480 


180 


17 


17.51(0.18) 


Feb 1995 - 


Oct 2000 


o 


9182 


260 


17 


17.45(0.11) 


Feb 1995 - 


Dec 2003 


P 


9739 


270 


23 


17.42(0.25) 


Feb 1995 - 


Feb 2000 



N is the number of images taken by the BATC telescope. 



In the right panel, we show the same field, as extracted 
from the BATC z-band images taken on 1995 Dec 20. 
SN 2004dj coincides with S96 to within 0.7" based on 
astrometry measurements in these two images. We ex- 




FlG. 1— (Left) An image of SN 2004dj in the BATC i band of 
the NAOC 60/90cm Schmidt telescope on 2004 Aug 11. The SN 
is circled. (Right) The same field is observed before SN 2004dj on 
1995 Dec 20. The SN is found to coincide with Sandage Star 96 or 
S96 (circled). The field of view of both images is 3' X 3'. 

tracted 213 images of NGC 2403, taken in 14 BATC fil- 
ters except for the b filter, from the BATC survey archive 
during Feb 1995 to Dec 2003. Table 1 contains the log 
of observations. Multiple images of the same filter were 
combined to improve the image quality assuming a non- 
variable S96. These serendipitous data provide a unique 
opportunity to examine pre-SN S96. 

2.2. Intermediate- Band Photometry of S96 

S96 is located in a somewhat complex background in- 
volving a spiral arm and possibly surrounding H II re- 
gions. To obtain proper photometry, we need to consider 
the flux contribution from the host galaxy background 
underneath S96. We invoked a method that assumes the 
spiral plane around S96 being stable and satisfying the 
Laplace equation and solves for the flux distribution of 
the spiral arm at the position of S96 (Zhang et al. 2004). 
The pure flux of S96 can then be estimated by subtract- 
ing the host galaxy contribution from the total flux. 

The final magnitudes of S96 are determined on the sub- 
tracted images with the standard aperture photometry. 
The BATC photometric system calibrates the magnitude 
zero level similar to the spectrophotomctric AB magni- 
tude system. For the flux calibration, the Oke-Gunn 



primary flux standard stars HD 19445, HD 84937, BD 
+26°2606, and BD +17°4708 were observed during pho- 
tometric nights (Yan et al. 2000). The results of well- 
calibrated photometry of S96 in 14 filters are summa- 
rized in the fifth column of Table 1 . These intermediate- 
band magnitudes (Fig. 3) agree well with the broad-band 
UBVI magnitudes of Larsen (1999). 

The numerous BATC images of NGC 2403 (Table 1) 
allow us to examine the variability of S96 in the past 
decade, which will shed light on the nature of the object. 
For example, significant light variations are expected on 
a time scale of a few years, if S96 were a luminous blue 
supergiant (Humphrey & Davidson 1994). We converted 
all BATC magnitudes, measured at different epochs and 
with higher qualities, into the i band by utilizing the 
color of S96 inferred from Table 1. As shown in Fig 2, 
the composite i-band light curve of S96 in the past decade 
does not show significant luminosity variations and most 
of the data are within ±0.13 mag (la) of the average 
value. The fit for a flat light curve gives x 2 = 49 for 43 
data points. This also justifies that the image combina- 
tion employed in our data reduction is reasonable. 
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Fig. 2. — The composite i-band light curve of S96. Solid circles 
show the i band magnitudes and other symbols are those converted 
from other 13 bands. All BATC band symbols are labeled on top. 
The solid line is the mean value of the i-band magnitudes and the 
two horizontal dashed lines mark ±1<t error. 

3. STELLAR POPULATION OF SANDAGE STAR 96 
3.1. Stellar Populations and Synthetic Photometry 

To probe the nature of S96, we compare its SED with 
theoretical stellar population synthesis models. We ex- 
plore two different SED families, one for single stars and 
the other for cluster populations. For stellar models, we 
used the theoretical stellar spectral flux library of Leje- 
une et al. (1997), covering a wide range of stellar param- 
eters. For cluster spectra, we used the SSP models (e.g., 
Bruzual & Chariot 2003; BC03), which assumed a single 
generation of coeval stars with fixed parameters such as 
age, metallicity and initial mass function. 

We convolve the SEDs of single stars from Lejeune et 
al. (1997) and stellar clusters from BC03 with the BATC 
filter transimision curves to obtain synthesized optical 
and near-infrared photometry for comparisons. The syn- 
thetic i-th BATC filter magnitude can be computed by 

m=-2.51og jA 48.60 1 

ipi(X)d\ 



Progenitor of SN 2004dj 
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where F\ is the SED from the library and ipi is the trans- 
mission function of the i-th filter of the B ATC photomet- 
ric system. Here, F\ varies with temperature and metal- 
licity for stellar models, and with age and metallicity for 
star clusters. We explore the nature of S96 by fitting the 
observed SED with distinct theoretical models. 

3.2. Reddening and Metallicity of S96 

To obtain the intrinsic SED of S96, the photometry 
must be de-reddened. The Na ID interstellar absorp- 
tion lines can provide clues to the reddening. However, 
the two reports about the measurements of these fea- 
ture lines of SN 2004dj (Patat et al. 2004; Gunther et 
al. 2004) show larger difference which prevents confident 
determination of the reddenings towards S96. As a re- 
sult, we treat E(B — V) as a fitting parameter determined 
simultaneously along with the effective temperature for 
stellar models or the age for cluster models. The values 
of extinction coefficient Ra for the BATC filters are ob- 
tained by interpolating the interstellar extinction curve 
of Cardelli et al. (1989). 

The SEDs for clusters or single stars are significantly 
affected by the adopted metallicity. Garnett et al. (1997) 
measured the metallicity and its radial gradient in NGC 
2403, using an expanded sample of H II regions. S96 is 
160" east and 10" north of the galactic nucleus of NGC 
2403, or at 4.0 kpc from the galactic center when re- 
moving the projection effect. At this radial distance, the 
relative oxygen abundance log(0/H) + 12 is 8.45 ± 0.10 
dex from Garnett et al., which is ~ 40% of the solar 
value. We then adopt a subsolar metallicity of z = 0.008 
for S96 in the analysis. 6 . 

3.3. Fitting Results 

We use the \ 2 test to examine which members of the 
SED families are most compatible with the observed one. 
We show in Fig. 3 the observed integrated photometry 
of S96 and the best-fit SED model. 

Our fitting results do not favor S96 as a single lumi- 
nous star. Fitting the BATC photometry for the stellar 
model yields the solution of a highly reddened \E(B — V) 
= 0.78±0.05 mag] B2IV star, but its best fit x is 24 for 
14 BATC points (with the DOF to be 12), indicating a 
rather poor quality of the fit. There are another three 
reasons why this fit does not work for S96: (1) the cor- 
rected intrinsic B — V color is too blue to be consistent 
with a B-type star; (2) such an extremely blue and bright 
star is most likely to be a luminous blue variable, but it 
did not show significant light variations (sec Fig. 2); (3) 
the high reddening would imply an unrealistically high 
luminosity, My < —18.9 mag at the maximum light, for 
SN 2004dj as a type II-P event (e.g., Hamuy 2003). 

In comparison, the star cluster models provide better 
fits to the observed SED of S96. The best reduced x 2 
of 0.76 is achieved for a cluster model with an age of 
19.1±3.3Myr and a reddening of E(B - y)=0.34±0.05 
mag. When the 2MASS infrared JHK S data (Skrutiskie 

6 We tested the validity of a presumed subsolar metallicity 
by treating it as a free parameter. The best-fit metallicity of 
0.012 ± 0.003 is obtained by using the BATC photometry and 
2MASS JHK S data (with a reduced x 2 ~ 0-72), which also prefers 
a subsolar value. However, uncertainties in the fitting for age and 
reddening increase significantly due to an additional parameter. 
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Fig. 3.— Upper panel: the best-fit SED (20 Myr solution) for 
S96, overlaid with that obtained by MA04 (dashed). The photo- 
metric points are shown with error bars (vertical ones for uncer- 
tainties and horizonal ones for the FWHM of BATC filter). Lower 
panel: predicted flux drops of S96 when the SN has faded away. 



TABLE 2 

Fitting Parameters for Sandage star 96 





BATC 


BATC+JHK S 


Age (Myr) 


19.1 ±3.3 


20.0 ±3.4 


E(B - V) (mag) 


0.34 ±0.05 


0.35 ±0.05 


Mass (10 3 M Q ) 


93.5 ±14.8 


96.2 ± 15.3 


Reduced \ 2 


0.76 


0.77 


Turnoff Mass (M ) 


12.2 


11.7 



et al. 1997) were included in the fit, the best-fit pa- 
rameters for S96 are a cluster age of 20.0±3.4Myr and 
a reddening of E(B — y)=0.35±0.05 mag with a simi- 
larly small reduced x 2 of 0.77. The reddenings yielded 
for these solutions are actually compatible with the one 
inferred from the reported Na ID feature lines (e.g., 
Patat et al. 2004), if the calibration by Munari & Z wit- 
ter (1997) is adopted. The uncertainties (l<r) in age 
and reddening are derived from the likelihood map that 
restricted by the x 2 distribution in the age-extinction 
plane. We reproduce in Fig. 3 the SED based on the 
parameters derived by MA04 (the dashed curve), which 
does not fit well the BATC data (especially in the i, j and 
k bands) since the best fit gives x 2 = 19 for 14 points. 

Adopting these age and extinction estimates, we de- 
termine the total mass of the stellar population from 
the distance to NGC 2403 and the measured V mag- 
nitudes (transformed from observed magnitudes in the 
BATC bands). We inferred the total mass of ~ 96,000 
M Q by comparing the measured V band luminosity with 
the theoretical mass-to- light ratios (BC03). The turn-off 
mass is ~12M Q for a cluster age of 20 Myr. The param- 
eters estimated for S96 are listed in Table 2. 

4. DISCUSSION AND CONCLUSIONS 

MA04 estimated a younger age of ~ 13.6Myr and a 
lower reddening 7 of E(B — V) = 0.17 mag for S96, us- 
ing the broad-band UBVI photometry of Larsen (1999). 

7 There is also a solution for an older age of 28.8Myr and a larger 
reddening of E(B — V) = 0.28 mag with a reduced x 2 of ~ 1-74 in 
their analysis but is considered less likely. 
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Our results differ from those of MA04 at a confidence 
level of ~ 2a, which is primarily caused by the differ- 
ences in metallicity and spectral coverage. Given a lower 
metallicity of z = 0.008 for S96, however, the combina- 
tion of UBVI and JHK S photometry yields two solu- 
tions in the age-reddening plane. The old solution of age 
~ 19 Myr with a reduced % 2 of 0.89 is similar to ours, 
while the young one of age ~10 Myr is equally compelling 
and cannot be rejected due to a smaller \ 2 of 1.02. The 
occurrence of multiple solutions may be related to the 
lower spectral resolution of the broadband photometry. 
In comparison, the BATC data constrain the SED better. 
For instance, our i-band photometry could indicate that 
S96 was not an H a -bright source and hence the young 
solution is unlikely, but the existing UBVI data alone 
fail to do that. The effect of different SED models on 
the results may be small as suggested by de Grijs et al. 
(2005). This is manifest in our data analysis of S96 us- 
ing both BC03 and Starburst99 models (Letherer et al. 
1999). At a fixed metallicity, the age difference derived 
from these two models is within 2 — 3 Myr and the ex- 
tinction difference is within 0.06 mag. 

While the observed SED of S96 resembles that of a star 
cluster, several photometric points shown in Fig. 3 do 
not fit well, e.g. the / and j bands. This small deviation 
of ~ 0.1 mag from the best-fit SED may be caused by 
the model uncertainties, the photometric errors, or even 
flux modulations of the immediate progenitor star. In 
the last case, we set limits on the possible variability. 
Given the luminosity of S96 and that predicted from the 
supergiant with a main-sequence mass of 12M Q (which is 
-6.9 < M bo i < - 7.8 mag; e.g., Smartt et al. 2003), we 
estimate that the flux contribution of the progenitor to be 
3 - 8% of the whole cluster, depending on the supergiant 
types. By relegating all scatters (~ 0.13 mag) shown in 
Fig. 2 to the progenitor itself, we place a rather crude 
upper limit of ~ 1.5 mag on its light variation. 
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The progenitor mass of SN 2004dj inferred in S96 
(~12M Q ) is within the mass range found for other type 
II-P SNe (Leonard et al. 2002; Van Dyk et al. 2003; 
Smartt et al. 2003, 2004). While it remains uncertain 
whether the progenitor of SN 2004dj was a blue, yellow, 
or red supergiant in the compact cluster S96, a compar- 
ison of flux changes for pre- and post-outburst phases 
of S96 in the blue and red bands may offer key clues to 
distinguish the various scenarios. Given a blue or yel- 
low supergiant, the post-SN flux of S96 would change 
by ~ 0.02 — 0.03 mag (Fig 3). If the progenitor was a 
red supergiant, the flux drop of S96 at NIR bands from 
k to p would reach ~ 0.07 mag. This change might be 
detectable by scrutinizing S96, preferably with the Hub- 
ble Space Telescope, when the supernova becomes signif- 
icantly dimmed several years later (Zhang et al. 2004). 

We conclude that SN 2004dj occurred in the young 
compact star cluster S96 of age ~ 20Myr. The lack of 
light variations in S96 supports this cluster identification. 
The progenitor of SN 2004dj is inferred to be a super- 
giant with a main-sequence mass of ~ 12M Q . Post-SN 
observations of S96 and detailed photometric and spec- 
troscopic studies of the SN evolution may further con- 
strain the progenitor of SN 2004dj. 
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